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I.  Introduction 


We  have  previously  demonstrated  the  usefulness  of  network  analyzers  for  plasma 
diagnostics  using  spheres  of  various  sizes.1'4  In  the  most  recent  work3’4  we  demonstrated 
a  method  of  determining  plasma  potential  and  the  electron  energy  distribution  from 
which  electron  temperature  can  be  estimated.  In  the  present  note  we  apply  the  same 
algorithms  to  a  cylindrical  probe  of  length  L  =  15  cm  and  radius  r  =  .16  cm.  This  work 
will  not  repeat  the  theoretical  framework  but  will  present  those  equations  used  in  the 
analysis  with  only  a  very  brief  description.  Interested  readers  are  referred  to  the  earlier 
works  which  form  the  theoretical  basis  of  the  technique  used.  '  ’ 


II.  A  Brief  Overview  of  the  Algorithm 


Electron  density  is  detennined  as  with  conventional  impedance  probes  by 
applying  a  small  rf  signal  and  detennining  plasma  resonance  in  the  real  part  of  the 
impedance  Re(Zac)  and  also  by  the  change  in  the  imaginary  part  Im(Zac)  from  capacitive 
to  inductive  ( i.e .,  a  zero  crossing  indicating  vanishing  reactance).1  As  an  approximation, 
the  bulk  plasma  behaves  as  inductance  whereas  the  sheath  region  near  the  probe  can  be 
characterized  as  an  ac  resistance,  Rac,  in  parallel  with  an  ac  capacitance,  Cac-  Since  the 
network  analyzer  provides  both  Im(Zac )  and  Re(Zac),  there  are  two  identifiers  of  plasma 
frequency. 


The  determination  of  the  electron  distribution  function,  /(e)34,  electron 
temperature',  Te  ,  and  plasma  potential  ,  <pp,  requires  that  we  operate  in  the  frequency 
range  copi  <  ca  <cope(r0),  where  oopi  is  the  ion  plasma  frequency,  co  is  the  applied  network 
analyzer  frequency,  and  cope(ro)  is  the  electron  plasma  frequency  at  the  surface  of  the 
probe  with  radius  ro.  Because  of  the  frequency  range,  the  ac  ion  current  is  small  and  the 
ac  resistance  is  given  by', 


d  t  dle  \-i 

(dvp} 


(1) 


where  Ie  is  the  dc  electron  current  at  bias  voltage  Vp  i.e.,  the  total  ac  current  is  now 
largely  electron  current.  The  Druyvesteyn7  equation  may  now  be  expressed  in  terms  of  a 
first  order  derivative  or, 


m= 


\me£ 

\dR;l  | 

J  2 

(2) 


and  thus  noise  is  reduced  in  the  calculation  of  the  distribution  function,  where  A„  is  the 
probe  area,  me  is  electron  mass,  and  s  is  the  energy.  In  addition,  earlier  work' 
demonstrates  the  relationship  between  the  real  part  of  the  analyzer  output  and  Rac, 
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(3) 


Re(ZflC)  = 
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1  +  (coRacCac) 


2  ‘ 


Therefore,  in  the  frequency  range  given,  and  provided  coRacCac  «  1,  we  have  Re(Zac)  ~ 
Rac  and  Re(Zac)  from  the  network  analyzer  output  can  be  directly  used  in  place  of  Rac  in 
Eqn  (2).  This  simplifies  finding  f(s)  by  avoiding  a  calculation  for  Rac  when  reducing  the 
data.  For  comparison  purposes,  we  use  both  methods  of  finding  f(s);  one  which  relies  on 
finding  Rac  and  the  other  which  uses  Re(Zac)  directly.  One ef(s)  is  known  we  can  estimate 
Te  assuming  a  Maxwellian  distribution  and  by  further  integration  find  another  estimate  of 
the  electron  density  for  consistency  comparison. 

Since  the  basis  for  determining  plasma  potential  is  that  d2Ie/dVp2  vanishes  at  Vp  = 


dRac 

=  -R2 

dVn 

ac 

dV; 

P 

vp=*p 

V  p  J 

and  finally, 

JRe(Zac) 


(4) 

(5) 


Plots  of  the  network  analyzer  output  of  Re(Zac )  versus  applied  bias  for  frequencies 
in  the  range  specified  will  then  show  a  minimum  at  Vp  =  (pp.  Due  to  the  dependence  of 
f(e)  on  dRac/dVp  seen  in  Eq  (2)  we  are  able  to  construct  f(s)  (Figure  3)  in  the  vicinity  of 
(pp  for  small  negative  bias  voltages  and  from  this  derive  ne  (Figure  4)  and  Te  (Figure  5)  as 
will  be  covered  below. 


III.  Outline  of  Experimental  Procedure 

We  refer  the  reader  interested  in  the  experimental  details  to  the  earlier  works  and 
only  provide  an  outline  of  that  same  description  here. 

The  experiments  were  conducted  using  as  a  probe  a  stainless  steel  cylinder  with 
length,  F  =  15  cm  and  radius,  r  =  0.16  cm  which  is  connected  to  an  HP8735D  Network 
Analyzer  through  50  Q  coaxial  cable  which  provides  the  driving  signal.  This 
arrangement  including  the  chamber,  analyzer  and  the  coupling  circuitry  is  shown 
schematically  in  Reference  2.  The  cylinder  is  mounted  on  a  1/4  inch  diameter  ceramic 
and  steel  support  which  is  connected  to  1/4  inch  diameter  semi-rigid  copper  50  Ohm 
coaxial  cable. 

The  determination  of  plasma  impedance  depends  upon  the  network  analyzer 
measurement  of  the  complex  reflection  coefficient,  r(co).  The  analyzer  returns  as 
separate  outputs  Re  Zac(co )  and  Im  Zac(co)  where, 


2 


z“w=z»[ri?T] 

1  -  ffiy) 


(6) 


and  Zo  (=50  Q)  is  the  internal  impedance  of  the  analyzer, 
reflected-to-total  power  is  given  by, 

2=^ 


We  also  note  that  the  ratio  of 


(7) 


where  P0  =  Pr  +  Pt  with  PR  and  PT  the  reflected  and  transmitted  powers,  respectively. 
(The  quantity  1-  \  r\  2  is  the  nonnalized  transmitted  power  and  this  output  is  also 
available). 

IV  Results  and  Discussion 


We  compare  two  data  sets  taken  under  similar  plasma  conditions  on  different  days 
with  the  probe  aligned  along  the  magnetic  field  direction.  The  magnetic  field  was  ~  2 
gauss  and  is  primarily  employed  to  contain  the  plasma  radially  by  inhibiting  diffusion  to 
the  walls. 

IV.  1  Electron  densities  from  Re(Zac)  and  Im(Zac) 

Figures  1  and  2  show  the  impedance  probe  results  for  the  2  cases.  The  Re(Zac) 
shows  resonance  at  the  plasma  frequency  and  the  Im(Zac)  shows  the  zero  crossing 
indicating  the  transfer  of  total  power  into  the  real  part  of  the  impedance  at  resonance. 
These  densities  are  compared  below  to  their  values  using  alternate  methods. 

IV.2  Magnetic  field  effects  in  cylindrical  Langmuir  probe  data  analysis 

With  the  probe  aligned  along  B,  magnetic  field  effects  under  certain  restrictions 
can  decrease  the  collected  electron  current  and  therefore  Langmuir  probe  determination 
of  density  based  on  electron  saturation  current  collection  indicates  an  artificially  lower 
bulk  plasma  value  than  is  the  actual  case.  Because  of  this,  electron  temperature 
determinations  can  also  be  compromised  since  most  fitting  routines  use  the  entire  IV 
characteristic  in  the  fit  itself.  The  effects  of  magnetic  fields  on  both  spheres  and 
cylinders  have  been  studied  over  the  years5’6,  and  it  is  not  the  intention  here  to  retrace  that 
work  although  some  general  comments  are  appropriate  when  comparing  to  impedance 
probe  measurements.  The  magnitude  of  the  effect  depends  on  the  electron  gyroradius,  re , 
in  addition  to  the  size  and  shape  of  the  probe  and  its  sheath.  In  a  thin  sheath  limit,  and 
assuming  the  same  applied  bias,  we  expect  generally  less  influence  as  a  function  of  probe 
orientation  with  respect  to  B  as  probe  radius,  r,  decreases  with  respect  to  re.  In  the 
present  case  the  sheath  radius,  rsh  is  estimated  at  0.5  cm  (assuming  rsh  ~  5 In)  while  the 
cylinder  radius  is  0.16  cm  or  rsh  >  r.  Also,  since  rjr  for  a  2  gauss  confining  field  is  of 
order  10  for  a  typical  Te  of  1  eV,  we  expect  less  of  an  effect  than  would  be  the  case  for  a 
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lower  ratio.  However,  it  has  been  shown  for  cylindrical  probes  that  even  when  re  >  r 
magnetic  fields  can  be  important  and,  in  addition,  even  extending  to  the  case  where  rsh  is 
comparable  to  re ,6  To  avoid  complications  associated  with  inferring  density  from  the 
electron  saturation  portion  of  the  characteristic  when  operating  in  a  magnetic  field  using 
conventional  analysis  techniques,  ion  saturation  current  is  often  used  instead  where, 


The  factor  e0'5  =  ee<p/Te  arises  from  the  density  decrease  at  the  sheath  edge  (where 
quasineutrality  applies  at  the  bulk  plasma)  to  the  pre-sheath.  With  the  ion  flow  velocity 
at  the  sheath  edge  taken  as  the  ion  sound  speed,  ion  energy  and  flux  conservation  requires 
approximately  that  e(p(rSh)  =  TJ2  and  the  Boltzmann  condition  then  allows, 

-ev^sk) 

nuJrsh)  =  >he  Te  =n 0W0'5  (9) 

Since  the  ion  gyroradius  is  typically  much  larger  than  probe  dimensions,  and  if  the  sheath 
is  small,  there  is  no  effect  on  Iion  from  a  large  range  of  magnetic  field  magnitudes.  The 
value  determined  however  may  still  be  compromised  by  the  structure  of  the  fitting 
routine  used  to  analyze  the  characteristic  i.e.,  density  detennined  from  even  the  ion 
saturation  level  can  be  incorrect  since  the  data  reduction  technique  typically  applies  a  fit 
to  the  entire  IV  characteristic,  a  portion  of  which  is  the  electron  saturation  current  region. 
Also,  the  existence  of  a  high  energy  electron  component  added  to  the  bulk  population  can 
complicate  this  analysis8  but  for  a  MB  distribution  the  analysis  is  justified. 


IV  3  -  Resonance  probe  determinations  of  Te ,  ne  using  f(s) 

In  addition  to  detennining  ne  with  the  resonant  probe  techniques  based  on  electron 
plasma  resonance,  we  are  also  able  to  detennine  plasma  potential,  Vp,  and  Te  based  on 
construction  of  the  electron  energy  distribution  as  outlined  above  and  in  cited  works.  In 
many  cases  of  practical  interest  we  are  simply  able  to  use  the  Re(Zac)  output  of  the 
network  analyzer  as  seen  in  Eqs.  (2)  and  (3)  when  the  approximation  that  ( (oRacCac )  «  1 
is  valid.  We  stress  that  Cac  seen  in  Eq.  (3)  is  detennined  from  experimental  values  with 
the  exception  of  an  assumption  of  a  constant  bulk  plasma  inductance  Lp4.  For 
comparison  purposes  we  show  in  each  plot  of  Figure  3  two  experimental  determinations 
of f(e)  using,  in  one  case,  Rac ,  and  in  the  other  Re(Zac).  These  cases  conespond  in  Eq  (2) 
to  the  situation  where  the  condition  above  is  not  rigorously  satisfied  and  in  fact  is  on  the 
order  of  1  itself  even  for  the  lower  frequency  ranges.  Citing  a  typical  value  we  allow  Cac 
~  2  x  10~u  F,f=  3  MHz  (oj=2nj)  and  Rac  ~  2000  Q.  In  this  example  (ojRmCacf  ~  0.56 
and  so  we  should  suspect  some  difference  in  the  construction  of  f(s)  depending  on 
whether  we  choose  Re(Zac)  or  Rac  in  Eq.  (2).  Examination  of  each  of  the  plots 
demonstrates  that  the  difference  is  minimal  and  only  begins  show  deviation  when  the 
applied  bias  becomes  larger.  As  a  result,  the  effect  on  the  detennination  of  ne  and  the 
average  energy  is  minimal  for  this  parameter  range  based  on  the  choice  of  method 
employed.  We  note  that  the  construction  of f(s)  depends  on  the  lrst  derivative  of  RacA  (or 
Re(Zac )~  )  and  not  on  the  function  itself. 
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Figure  4  shows  a  comparison  among  the  three  methods  of  determining  ne  from: 
(1)  the  Langmuir  probe  sweep,  (2)  the  resonance  peak  and,  (3)  integration  of f(e).  Figure 
5  shows  comparison  of  the  Langmuir  probe  sweep  to  detennine  Te  versus  estimating  it  as 
Te  =  2/3  Eav  where  Eav  ,  the  average  energy,  is  found  buy  integration  over  ef(e).  (Note 
that  the  abscissa  of  Figures  4  and  5  is  applied  frequency.  The  reason  for  this  is  for 
consistent  representation  of  measurements  taken  using  the  Langmuir  probe  sweep,  i.e.,  a 
value  for  ne  ,  Te  was  determined  before  each  sweep  of  the  bias  voltage  for  a  given 
frequency.  Since  the  other  methods  determine  a  single  value  only,  they  appear  as  a 
straight  line  on  each  plot.)  Comparing  both  cases  of  the  resonance  peak  determination  of 
ne  shows  a  higher  value  as  would  be  expected  from  the  Langmuir  probe  analysis  if  there 
is  a  magnetic  field  effect  on  sheath  collection  due  to  orientation  as  discussed  above, 
However,  integration  of  the  distribution  function  yields  a  value  closer  to  that  obtained  by 
the  Langmuir  probe  sweep  analysis.  The  reason  for  this  is  not  immediately  evident.  One 
possible  explanation  lies  in  the  fact  that  the  construction  of  f(e)  is  limited  insofar  as  the 
range  of  bias  voltage  covered  due  to  experimental  noise  limitations  at  lower  frequencies. 
This  suggests  we  are  “missing”  a  tail  of  the  distribution  further  suggesting  that  the 
integration  of f(e)  yields  a  value  less  than  is  the  actual  case. 

IV.  4  Re(Zac )  vs  Vbias  at  fixed  frequencies 

Finally,  Figure  6  shows  Re(Zac)  vs  Vbias  at  fixed  frequencies.  As  described  briefly 
above  and  in  the  works  cited,  this  is  the  primary  plot  from  which  plasma  potential  is 
obtained  and  from  which  the  values  of  Re(Zac)  arise  for  use  in  the  construction  of f(e). 
The  local  minima  at  approximately  2  volts  in  this  figure  are  consistent  with  the 
determination  of  Vp  from  the  Langmuir  probe-based  characteristic. 


V  Summary 

We  have  developed  an  impedance  probe  technique  for  determining  electron 
density,  temperature,  plasma  potential  and  the  electron  distribution  function  over  the  past 
three  years  in  studies  using  as  rf  impedance  probes  spheres  of  different  sizes.  The 
theoretical  basis  for  this  work  found  in  the  references  in  Section  VII  suggests  that  the 
technique  is,  to  a  large  degree,  independent  of  probe  geometry.  The  work  presented  here, 
using  a  cylindrical  probe  with  a  high-aspect  ratio,  is  consistent  with  this  conclusion. 
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VI.  Figure  Captions 


Figure  1  -  The  real  part  of  the  complex  plasma  impedance,  Re(Zac),  plotted  versus  the 
applied  network  analyzer  rf  frequency  for  two  separate  experimental  series.  The  primary 
resonance  occurs  at  the  bulk  plasma  electron  frequency.  Resonances  of  lower  amplitude 
are  produced  by  the  chamber  wall  reflections  and  are  unrelated  to  the  bulk  plasma. 

Figure  2  -  The  imaginary  part  of  the  complex  plasma  impedance.  Im(Zac),  plotted  versus 
the  applied  network  analyzer  rf  frequency  for  two  separate  experimental  series  using  a  1 5 
cm  x  0.16  cm  (length  x  radius)  cylinder.  The  zero  crossings  occur  at  the  bulk  plasma 
electron  frequency.  As  with  the  Re(Zac)  structures  in  either  the  capacitve  or  inductive 
reactance  are  associated  with  chamber  wall  reflections  and  are  unrelated  to  the  bulk 
plasma. 

Figure  3  -  Plots  of  the  electron  distribution  function  f(e)  vs  applied  probe  bias  for  an 
analyzer  sweep  frequency  of  3  MHz  as  determined  by  2  separate  methods:  (1)  The  black 
curve  shows  f(e)  using  Re(Zac)  in  Eq.  (2)  and,  (2)  the  red  curve  shows  f(e)  using  Rac  in  Eq. 
2.  The  significance  of  this  is  described  in  the  body  of  the  text. 


Figure  4  -  The  bulk  electron  plasma  density,  ne  ,  as  determined  by  3  separate  methods 
indicated  on  the  figure:  (1)  The  data  points  are  for  separate  voltage  sweeps  using  the 
cylinder  as  a  Langmuir  probe,  (2)  The  red  line  is  the  value  of  ne  as  detennined  by  the 
resonances  shown  in  Figures  1  and  2  and,  (3)  The  green  line  is  the  value  of  ne  as 
determined  by  integration  of  the  distribution  function  shown  in  Figure  6.  Construction  of 
the  distribution  function  is  described  in  the  body  of  the  paper. 


Figure  5  -  Plots  of  the  electron  temperature,  Te  ,  as  determined  by  2  separate  methods 
indicated  on  the  figure:  (1)  The  data  points  are  Te  as  determined  by  conventional 
Langmuir  probe  analysis  and,  (2)  The  red  line  is  Te  as  determined  by  integration  of  the 
electron  distribution  function  to  obtain  the  average  energy.  Construction  of  the 
distribution  function  is  described  in  the  body  of  the  paper. 

Figure  6  -  Plots  of  the  Re(Zac)  for  frequency  scans  varying  from  2  to  9  MHz  vs  applied 
bias  voltage  for  two  separate  experimental  series.  For  both  runs  the  plasma  potential 
determined  from  Langmuir  probe  data,  (pp~  2  V,  corresponding  to  the  local  minimum  at 
this  bias  in  both  plots  shown. 
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Ne  vs  frequency  - 100-1  cylinder  -  022510 
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Figure  4 
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Te  vs  frequency  -  1 00-1  cylinder  -  022510 
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T  vs  frequency  - 100-1  cylinder  -  030210 
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Figure  5 
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ReZ  vs  Vbias  - 100  to  1  ratio  cylinder  -  022509  run  1 
(f=  2  to  9  MHz) 
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